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We show how correlation functions of the spin- 1/2 Heisenberg chain without magnetic 
field in the anti-ferromagnetic ground state can be explicitly calculated using information con- 
tained in the quantum Knizhnik-Zamolodchikov equation [qKZ]. We find several fundamental 
■ relations which the inhomogeneous correlations should fulfill. On the other hand, it turns out 

that these relations can fix the form of the correlations uniquely. Actually, applying this idea, 
we have obtained all the correlation functions on five sites. Particularly by taking the homoge- 
neous limit, we have got the analytic form of the fourth-neighbor pair correlator (5|S , | +4 ). 
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1 Introduction 



The anti-ferromagnetic spin- 1/2 Heisenberg XXX chain d 

L 

H = (Sj Sj +1 + SjSj +1 + SjSj + ^j , (1-1) 
i=i 

is one of the most fundamental models in the study of quantum magnetism in low dimensions (see, 
for instance, the book |2J). 

The Hamiltonian (11.11) was diagonalized by Bethe in 1931 [3| and the ground state in the 
thermodynamic limit was investigated by Hulthen in 1938 H. The physical content was clarified 
by Faddeev and Takhtajan @ using the algebraic Bethe ansatz formulated by Faddeev, Sklyanin 
and Takhtajan (see the review J5J). Exact integrability of the Heisenberg model is intimately 
connected with the Yang-Baxter equation 0. 

The Heisenberg model also appears as an effective Hamiltonian of more complicated systems 
such as the one-dimensional Hubbard model at half-filling with a large coupling constant U. 

Since discovering of the model a lot of physical quantities have been calculated exactly which 
afford significant and reliable data for the experiments as well as numerical simulations. However, 
most of them are the bulk quantities and we have obtained little results on the correlation func- 
tions. Especially significant are the (static) spin-spin pair correlators (SjSj +k ) (or equivalently 
(S^ Sj +k ) /2; Sj = Sf ± iSj), for which only the first and the second neighbor (k = 1, 2) have 
been known for a long time: 

(S?S? +1 ) = — - -In 2 ~ -0.14771572685, (1.2) 
J J 12 3 

(S*S* +2 ) = ^-Un2 + |c(3) - 0.06067976996. (1.3) 

Here £(s) is the Riemann zeta function and (■ ■ •) denotes the ground state expectation value. Note 
that (11.21) was derived directly from the ground state energy in the thermodynamic limit by Hulthen, 
while (11.31) was obtained by Takahashi in 1977 via the strong coupling expansion for the ground 
state energy of the half-filled Hubbard model [ 8 1 (see also another derivation by Dittrich and In- 
ozemtsev [9|). Unfortunately we cannot generalize these methods so as to calculate further corre- 
lations (S?S? +k ) k > 3 . 
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On the other hand, in 1990s, correlation functions were studied systematically for the Heisen- 
berg XX Z model. It was mainly developed by Kyoto group (Jimbo, Miwa, Miki, Nakayashiki, 
lfi"4l . see also the book [15| ). Utilizing the representation theory of the quantum affine algebra 
U q (sl2), they first derived multiple integral representation of arbitrary correlators for for the mas- 
sive XX Z anti-ferromagnet. In 1996 Jimbo and Miwa extended the above results to the XXX 
chain and the massless XX Z chain [17]. The manifest integral formula for special case of corre- 
lation functions, Emptiness Formation Probability (EFP) of the XXX model was earlier obtained 
by Korepin, Izergin, Essler and Uglov in [ 16 1. Strictly speaking the formulas in the massless region 
were a conjecture since they were obtained as solutions to the quantum Knizhnik-Zamolodchikov 
equations (qKZ equations) with a certain assumption. However, the same integral formulas have 
been reproduced by Kitanine, Maillet and Terras, ll20ll2T1 in the framework of the Quantum Inverse 
Scattering Method. So we have now at hand the multiple integral representations for any arbitrary 
correlation functions for the Heisenberg XX Z chain. 

In 2001, Boos and Korepin invented a method to evaluate these multiple integrals for the 
Heisenberg XXX chain ll22ll23ll . They have shown that the integrand can be reduced to a certain 
canonical form so that the integrations can be performed. They applied the method to the case of 
the Emptiness Formation Probability. There has been also suggested conjecture that the correlation 
functions for the XXX model may be expressed in terms of the Riemann's zeta function at odd 
arguments with rational coefficients. For instance, the results for P(4) and -P(5) look as follows: 

P(4) = ((1/2 + (1/2 + S] +1 ) (1/2 + fi» +a ) (1/2 + 5| +3 )> 

= I _ 2 In 2 + — C(3) - — In 2 ■ C(3) - -C(3) 2 - -C(5) + - In 2 ■ C(5) 
5 60 SV ; 6 SV ; 80 SV ; 24 sv ; 24 sv ; 

(1.4) 

P(5) = ((1/2 + 5-) (1/2 + S* +1 ) (1/2 + S» +2 ) (1/2 + S; +3 ) (1/2 + 5J +4 )) 

1 10, 281 > . , 45, A/ , 489 ^ x2 6775^ , 1225, 

= 6 " T ln 2 + !T C(3) " T ln 2 ' m ~ l6 Ci3) ~ T5J C(5) + 1" ln 2 ' C(5) 

_ ^ c(3 ) . c( 5) - HH5 C(5)2 + 5!5 C( 7) _ 1^5 ta 2 • «7) + ^CO) ■ cm 

g 4 SV i SV ! SV ) 25g SV ) g 4 SV I SV ! SV J 

(1.5) 

The last result was obtained in [28 1 . 

More recently the method was applied to calculate other correlation functions of the Heisenberg 
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chain dl.lt . Here let us list the results in [ 29 1 where all the independent correlation functions on 

four lattice sites are obtained: 

1 37 14 3 195 95 

<^ +3 ) = - - 3 In 2 + ^C(3) - y In 2 ■ C(3) - §C(3) 2 - ^C(5) + y In 2 ■ C(5) 

= -0.05024862725- ■• , 

I 1 Q1 1 3 35 5 
(S]S* +1 S° +2 S> +3 ) = - + T2 in 2 - - C (3) + - in 2 . ((3) + _ { (8)> + -«5) - - fa.2 ■ «5) 

II 77 5 3 65 5 

(SfS^Sf^S^) = In 2 + — C(3) - — In 2 • C(3) - — C(3) 2 - — C(5) + - In 2 • C(5) 

\ 3 j+i j+2 j+3/ 240 6 120 12 w 20 96 6 

11 169 5 3 65 5 

(S*S*S*S*) = — - - In 2 + — C(3) - — In 2 • C(3) - — C(3) 2 - — C(5) + - In 2 • C(5) 
\ 3 j+i j+2 j+3/ 240 4 240 12 w 20 96 6 

(1.6) 

Note that the other correlation functions among four lattice sites are expressed as a linear combi- 
nation of the correlations above. For example we have 

/ QZ QZ QZ QZ \ / QX QX QZ QZ \ _|_ / QX QZ QX QZ \ I / QX QZ QZ QX \ 

^+1^+2^+3/ - \ D i ^+1^+2^+3/ + \ D j ^'+1^+2^+3/ + Wi ^+1^+2^+3/ 
1 1 99 9 91 95 35 

= _L _ L l n 2 + -((3) - - In 2 • C(3) - -C(3) 2 - — C(5) + — In 2 • C(5), 
80 3 30 SV ; 3 SV ; 80 SV ' 96 SV ; 24 sv ; ' 

13 1 1 

P (4) = — + - (SjSj +1 ) + - (SjSj +2 ) + - (SjSj +3 ) + (SjSj +1 Sj +2 Sj +3 ) , 

(1.7) 

etc.... 

In papers ll30l I3T1 l32ll the above results were generalized to the XX Z Heisenberg chain both 
in massless and in massive regime. 

Generalizing the results (1 1 .6b to the correlation functions on five lattice sites is a natural next 
problem. In principle, it is possible to calculate them similarly from the multiple integrals as was 
actually done for P(5). It, however, will take a tremendous amount of work to complete them. 

It is worthy to note that all the above results have one specific feature, namely, they all are 
expressed in terms of one-dimensional integrals. In series of papers Il25ll26ll27l this phenomenon 
was explained by means of the above mentioned connection of correlation functions in the inhomo- 
geneous case with the quantum Knizhnik-Zamolodchikov equation, more precisely, by means of a 
duality between solutions to the qKZ equation on level -4 and 0. It was also argued that this fact 
has a deep mathematical origin based on theory of deformed hyper-elliptic integrals, symplectic 
group and special type of cohomologies investigated earlier by Nakayashiki and Smirnov lfT8lfT9ll . 



One of basic outcomes of the above scheme is a general ansatz for correlation functions in which 
only one transcendental function participates with some rational functions as coefficients. Some 
special additional properties derived from the qKZ equation allowed to fix P{n) until n = 6. 

In recent paper ll33ll the properties of the qKZ equation have been further investigated and the 
above ansatz has been proved by means of some special recursion relations which linearly connect 
correlation functions at n, n — 1 and n — 2. Coefficients in these relations appeared to be related 
to special transfer matrices over an auxiliary space of "fractional dimension". The solution to 
these relations has also been found in the paper ll33l which gives an explicit result for rational 
coefficients before the transcendental functions staying in the ansatz for correlation functions. In 
spite of the fact that this formula is explicit for any number of sites n, it is still hard to use it for 
getting manifest results for concrete n. 

In this paper we generalize the scheme suggested in paper [25| to the case of an arbitrary 
correlation function. We use a set of relations for correlation functions that follow from the qKZ 
equation. We optionally call them "first principle" relations. The claim is that these relations 
together with the above mentioned ansatz completely fix rational functions which participate in this 
ansatz. We find the whole set of correlation functions in case of n = 5. After taking homogeneous 
limit we come to the result in terms of the Riemann zeta function in full agreement with the above 
conjecture, which now should be considered as proved. The "first principle" relations as well as 
results for correlation functions should be definitely equivalent to those given by the formula (3.22) 
in the paper [33]. But it is still to show such an equivalence explicitly. 

2 Quantum Knizhnik-Zamolodchikov equations and fundamen- 
tal relations for correlation functions 

In this section we would like to describe an algebraic scheme which allows to calculate any corre- 
lation function for the XXX model. In fact we generalize here the method applied in paper lf24l 
to the particular case of the EFP. 
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Let R(l) be the i?-matrix for the XXX model 
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(2.1) 



(2.2) 



where it is implied that I = nz. 

The corner stone of our scheme is the quantum Knizhnik-Zamolodchikov equation [qKZ] ifTOl 
[T31 ITT1 [T2l . It is a linear system of equations for a vector function g n . In particular, the qKZ 
equation on level —4 looks as follows 



R(lj — lj + i)J'j +1 g n (h 



, /j, lj+l, , ^2n)ei, -- , 



(2.3) 



9n{h,--- ,kn-l,hn + 2lTi) tu ... tt2n _ ltt2n — g n (hn,h, - • ' , ^2n-l)e 2 „,ei,- ,e 2 „_i (2-4) 

It was suggested by the Kyoto school ifbUl that the correlation functions are connected with solu- 
tions to the qKZ equation. In particular, for the massless case Jimbo and Miwa ifTTl found that one 
should take a special solution of the qKZ on level -4 which satisfies one additional constraint 

= tfitj 9n-l(h, ■ ■ • , lj+2, ■ ■ 4 , ^n)e ll -,ej-i,e j+ 2 1 -,£ 2n (2-5) 
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Solutions to the above equations (12.412.51) are meromorphic functions with possible singularities at 
the points 

Sfo - l k ) =ttI, I e Z\0 

but regular at $s(lj — l k ) = ±ir. 

Due to Jimbo and Miwa, any correlation function 

pe' lr ..,t' n _ / eVj TTKE^X 
ei,...,e n \ ei ' " " e„ I 

is given by the following formula 

P et--i ■■-,Z n )= 9n(h, ...,l n Jn + ™, h + 7™) ei| ..., £n _ e / n) ... (2-6) 

where as above lj = nzj . 

Let us list the "first principle" relations for the correlation functions which follow from the 
qKZ equation: 

• Translational invariance 

Peti (z 1 + x,...,z n + x) = P<;;;;;f: (z u ...,z n ) (2.7) 

• Transposition, Negating and Reverse order relations 



Pei,...,e n . . • , Z n ) — P , ' ' , ( Z\, . . . , Z n ) — 

P-ell..'-e^{ Z li ■■■iZ n ) — PeZ-'ei( — Z n , . . . , — Zl) (2.8) 

• Intertwining relation 



K e'e' 

3 3 

P ' 



7 + 1 

ej,gj+i.. 



3; 



j+1 



IP. 



■ e 3 + li e J-- 



. z 



Jl 



Zj+i ■ ■ ■)RLl + + r {zj - Zj+i) 



(2.9) 
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• First recurrent relation 



> c 2> 
£2,- 



K z + i, z, z 3 . . . , z„ 



A c' c p e 2> e 3>" 

■^1 ,-e2 e i e 2-^_ e / i)e3i 



-Pel, £2, ■■■,&n h Z i %3 



K z, z 3 . . 
fa z 3 . . 



• 1 %r, 



(2.10) 



• Second recurrent relation 



lim Pel].. 



fa, 



-3 1 ■ ■ ■ 1 



Zn) = S ej;€ r ~ Pi} 



fa, 



-31 ■ ■ ■ 1 



Zn) (2.11) 



• Identity relation 



fir 



'£(zi, ■ ■ ■ ,z n ) 



Cl, • • • , £n 

Si e i = Si € i 



1 c ra 



Si e i = Si e « 



= p+;:::;+fa, ■ ■ ■ , z n ) = p~;:::;~ fa, ...,z n ) (2.12) 

Reduction relation 



p-t-> e 2: 



■ "fa, z 2, 



-l,+ fai Z 1i 



+ 



p ,t 2 ,...,t n 

r — ,e2,—,e n 



pz: 



n-V 



2l, ^2; • • 

-fa, Z 2, 



pe 2 ,. 

ea,.. 



r ei,-. 



(2.13) 



The relations (12.71) . (12.81) and (12.131) were discussed by Jimbo and Miwa in the paper ifTTll . The 
intertwining relation (12.91 ) is the direct corollary of the qKZ equations. The proof of the first and 
the second recurrent relations d2.10l2.lll ) can be performed in the same lines as it was described 
in the paper [24|. In particular, the relation (12.1 II ) was completely proved in the paper [33 1. The 
identity relation (12.121) is a direct consequence of the sl 2 symmetry of the solution to the qKZ. 
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In papers Il2"51l2"6l . a general ansatz for the correlation functions was suggested. Then it was 
proved 1 in ff33l 

[f] 

Yl A i]:::,t( k i, ■ ■ ■ ,k2m\zi,...,z n )G(z kl -z k2 )---G{z k2m _ 1 -z k2m ) 

(2.14) 



where 

CO 

G(*)= 2 £(-!)*• Jfc.£_nA (2.15) 



x 2 + 1 

x* + k 2 
k=i 



and At{]Z',eZ(ki, • • * , ^7m\ z it ■■■■> z n) are rational functions with known denominator and polyno- 
mial numerator of a known degree in variables z\,...,z n , namely, 

Ae' ...V n i I \ Qei,...,e n {klt ' ' ' , fem \ z l, ■ ■ ■ , Z n ) « ln 

A e] 'e:(h, ■ • • , feam|2l, • • • , «n) = / v ( 2 -16) 



ei,...,e 



where n^j^ ~~ z i) means tne product in which the differences with {i,j} = {k 2 i, k 2 i + i} and 
i,j(£K should be excluded, K is the set ki, - ■ ■ , k 2m . QtuZ',ln(ki, ■ • • , k 2m \zi, . . . , z n ) are poly- 
nomials of the variables z which have the same degree in each variable as in the denominator and 
also the same degree of homogeneity. Some useful properties of the function G(x) are listed below 

G(-x) = G(x), (2.17) 

G(±oo) = -~ (2.18) 

G(x-i) = a(x) +j(x)G(x), (2.19) 

where 



x — 2i x(x — 2i) 

a(x) = r , -){x) = \ ' (2.20) 

x — l x z + 1 



Moreover G(x) is a generating function of the alternating zeta functions 

CO 

G(x) = -2(1 + x 2 ) ]T(-1) k x 2k ( a (2k + 1), (2.21) 



k=0 



'The function cu(l) used in 1 33 1 is related to the function G by = G(il) + 1/2 



with 

( a (s) = = (1 - ^~ S )C(s). (2.22) 

n=l 

Conjectural formula for the very first term when m = looks as follows 

Atzii^i, ...,z n )= Qi';;±{^ ■ ■ ■ . *0 = (2-23) 

where r is the number of down (or up) spins in the set ei, . . . , e n (or e[,..., e' n ). 

Our claim is that the "first principle" relations (I2.7I) - (I2.13I) together with the ansatz (12.141) and 
conjecture (12.231) completely defines all unknown coefficients in polynomials 
Qel','.'.'.feZ(ki, • • • , k 2 m\zi, ■ ■ ■ , z n ) and hence the whole set of correlation functions. 

Here let us list the explicit form of the inhomogeneous correlaion functions for n = 1 and 
n = 2. 

P+(z l )=P = (z 1 )= 1 -, 

P++{zi, to) = P~(z u z 2 ) = - + -G{z x - z 2 ), 

p+:( Zl , z 2 ) = p:+( Zi , z 2 ) = \- l -G{ Zl - Z2 ), 

DO 

P-+( Z1 , z 2 ) = P++( Zl , z 2 ) = l - + ^G( Zl - z 2 ). (2.24) 

One can easily confirm these correlation functions satisfy the first principle relations. 
In the subsequent sections, we denote the inhomogeneous EFP as 

Pn( Z\-, z 2 , z n ) = P^tf.'" ) +(^i, z 2 , z n ) (— P (zi, z 2 , z n )) . (2.25) 

Especially we have P 2 (z 1 , z 2 ) = P++(zi, z 2 ). The explicit form of the above inhomogeneous EFPs 
up to n = 6 was obtained in ll2"4ll . 

3 Inhomogeneous correlation functions for n = 3 

Here we consider the inhomogeneous correlation functions for n = 3 and the corresponding first 
principle relations. This case was considered as an example in the paper ll33l . But we reproduce 
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p e 'l e 2 £ 3 


- + + 


+ - + 


+ + - 


- + + 


To H~ TniGo^ — G-\o — Cms) 

12 12 v ^° t 
+ 12^123 


12 + 6^12 + 1^123 


12 10 

i riA 
12^123 


+ - + 


1 _i_ 1 ^ — * i f~iA 
12 G^ 12 12^123 


To 4" To (G13 — G*23 — G12) 
-TJC123 


12 + 1^23 + 7^123 


+ + - 


12 + \ G \3 - |Ci23 

1 i r<A 
12^123 


1 _i_ 1 / — * i r<A 
12 "T 6^23 12^123 


12 + T2~(^ 12 — ^13 — C23) 
+ nG*i23 



Table 1: Inhomogeneous correlation functions forn = 3 in the sector t\ + e 2 + e 3 = 1 



it here in order to make our approach more illustrative. To describe the explicit expressions of 

e f e' e f 

P £1 1 £2 2 (E3 3 (2; 1 , z 2 , z 3 ), we introduce some notations. 



Gi 



r*iA 
^123 



n _ G\2 G13 G 23 
^123 = 1 I 



G(Zij), 

z 13 z 23 Z 12 Z 32 z 2\ z 3\ 

— -—) Gi 2 +(- -)g 13 +(- j G 2 ,. 

Z XZ z 23/ \ z 23 Z Y1 J \ z 2,\ z 21, 



(3.1) 



where we have also introduced an abbreviation z^ = Zi — Zj. With these notations, the EFP 

P 3 (z u z 2 ,z 3 ) = P^(z u z 2 ,z 3 ) is given by J24| 



-^3(^1, z 2, z 3) — J + T7j ^ 13 



C23) — Y^Ci23- 



(3.2) 



Next we consider the inhomogeneous correlation functions Pei*£ei{zi, z 2 , z 3 ) in the sector 



ei + e 2 + e 3 



e' + e' 



1. For convenience, we first exhibit the final results in Table 1. 



There each row corresponds to the subscript e 1 , e 2 , 63, while each column corresponds to the su- 



perscript e'^e^eg. 



Below we show the strategy to calculate these correlation functions from our first principle re- 
lations. Assume we have already obtained the EFP, P 3 (zi, z 2 , z 3 ) and all the correlations for n = 2. 
As is shown in Table 1, there are 9 correlation functions in the sector t\ + e 2 + e 3 = e[ + e' 2 + t' 3 = 1. 
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However, we find, they are actually connected through the intertwining relation ( 12.91 ) as follows, 



P+-+{zi, Z 2 , z 3 ) — 2 3) 2 2 ) Z 2 , Z 3 ), 

z 23 Z23 

P++-{%U z 1i Z3) — P++- ( Z 2, ^i, Z 3 ) P + +l(Zi, Z 2 , Z 3 ), 

Z21 Z21 

P-++{ z li z 2, z 'i) = P+-+ ( z 2, Z\, Z 3 ) H {P_++ Zi, Z 3 ) — P_++(Zi, Z 2 , 2 3 )} , 

Z21 

P-++ Z 2i Z3) = — P-++{Zl, Z 3 , Z 2 ) P-++ ( Zlj Z2j Z3 )i 

z 32 z 32 

P+-+ (zi, z 2 , z 3 ) = P_++ (z 2 , Zi, z 3 ) H {P+-+ ( z 2, z\, z 3 ) — P + -+ (zi, z 2 , z 3 )} , 

Z21 

P+-+ ( Z l, z 2, Z 3 ) — P-++ ( Z 2, Zi, Z 3 ) P-++ { Z l> Z 1l Z z)l 

Z\2 z \2 

P++r(zi, z 2 , z 3 ) = P+I+( z ^ *a> z *) + — {K+- ( z u *3, z 2) ~ Pt^{z u Z 2 , Z 3 )} . (3.3) 

Z 3 2 

These relations imply we need to calculate only two correlations, for example, P~++(zi, z 2 , z 3 ) 
and P++*(zi, z 2 , z 3 ). The remaining correlation functions are obtained immediately from ( 13.31) . 
Moreover the diagonal one, Pr++ z 2 , z 3 ) can be derived easily from the reduction relation 

dm, 

PZXt{zu z 2 , z 3 ) = P 2 (z 2 , z 3 ) - P 3 {zi, z 2 , z 3 ) 

= g + g^*23 — + ^ (^ 12 + ^13 + G 23 ) — Y2^ 123 

= h + h (G23 " Gu ~ Gis) + h G ™- (3 - 4) 

Therefore our task reduces to calculate one off-diagonal correlation P++^(zi, z 2 , z 3 ). Note that 
there are different ways to connect the correlations from (13 .3b . They, however, are consistent with 
each other due to the Yang-Baxter relation for Rij(z). 

Look at other first principle relations and apply them to P++*(zi, z 2 ,z 3 ). First, we find the 
reverse order relation gives 

P-++^ 1; z 2 , z 3 ) = P+++(z 3 , z 2 , Zl ), (3.5) 

which means P++*(zi, z 2 , z 3 ) is symmetric with respect to the exchange z\ <-> z 3 . Next, the trans- 
position relation P^^(z 1 , z 2 , z 3 ) = P++l{—Zi,—z 2 ,—z 3 ) together with (13.31) . (13.41) and (13.51) 
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gives an equation 



= i z U z 2-, z 3 ) P++-( z 2, z li Z 3 ) P + i-{zi, Z 3 , Z 2 ) 

^12 ^32 #12 #31 £32 ^13 

' -P 3 (*i,*i) + — ^ ± ip a (z 1 ,^) + — ^±^2(^3)- (3.6) 



^12^32 z 31 z 32 Z13 Z\ 2 

The recurrent relations for Pl_+I(zi, z 2 , z 3 ) are given by 

P+t-( z + h z, z 3 ) = 0, P-++(z - i, z, z 3 ) = -P;+(z, z 3 ), (3.7) 
lim P7±±(z 1 ,z 2 ,z 3 ) = lim P^(z 1 , z 2 , z 3 ) = 0, 

21— >oo 23 — >oo 

lim P+ 22, z 3 ) = Jp+- 23). (3-8) 
The identity relation gives a relation, 

^±lp-++(z 1 , Z2 , z 3 ) + ^i±l P -++( 2l) Z3> Z2 ) = 2 P 3 ( Zl , z 2 , z 3 ) - P 2 (z 2 , z 3 ). (3.9) 

^32 ^23 

where we have used the relations (13.31) and (13.41) . Equations (I3.5M3.9I ) with the translational 
invariance 

P-++(2 1 + x, z 2 + x, z 3 + x)= P7+Z{ Zl , z 2 , z 3 ). (3.10) 

are our first principle relations reduced to P+?t(zi, z 2 , z 3 ). 

According to the ansatz (12.141) and (12.231) . let us assume P++7l(zi, z 2 , z 3 ) in the form 

TJ-++, \ 1 Cn + C\ 2 Zi 3 + Ci 3 z 23 + CuZi 3 z 23 

[Zl, Z 2 , Z 3 ) - — H — — Cri2 

Z\ 3 Z 23 

C21 + C22^i2 + C23^32 + C24^12^32 n 

H ^13 

Z\%Z 32 

. c 31 + C32^21 + C33Z31 + C 3 4Z 2 iZ 3 i 

H ^23, (J- 11 ) 

^21^31 

where are the coefficients which will be determined by the relations (I3.5I) - (I3.9I) . Note that from 
the symmetry relation (13.51) . we at once have several relations 

C31 = Cn, C32 = C13, C33 = C12, C34 = C14, C23 = C22. (3.12) 
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Relations (I3.6I H I3.9I ) generate equations for the coefficients c^-. Although they are an overdeter- 
mined system of the equations for c^, one can figure out it has a unique solution 

cn = c 2 i = C31 = C12 = -C13 = c 22 = c 23 = -C32 = c 33 = -J2, 

cu = c 34 = 0, c 24 = |- (3.13) 

Namely the correlation function P^j^i, 2 2 , ^3) is determined as 

6 [ 212232 2 V232 212 / J 



+ ids - + X (3.14) 



6 [ ^21^31 2 \Z 3 i 221/ 

= -)- — Gri3 — — Cri93 ~i~ — 

12 6 6 12 

Other correlation functions are obtained from the intertwining relations (13.31) . In this way we 
could calculate all the inhomogeneous correlation functions as shown in Table 1. The correla- 
tion functions in the sector e\ + e 2 + e 3 = e' x + e' 2 + e 3 = — 1 are simply obtained by the negation 
relations, 

piti&{z x , z 2 , Z3 ) = p::i;::t::'!(z u z 2 , z 3 ). (3.i5) 

Now let us think of an inhomogeneous generalization of the next-nearest neighbor correlation 
function, 

(S'SI) {z u z 2 , z 3 ) = - {P 3 {z u z 2 , z 3 ) - P++r(zi, z 2 , z 3 ) + P+r|(zi, z 2 , z 3 ) - Pr++(zi, z 2 , z 3 )} 
= i - P 2 {z 1 , z 2 ) - P 2 (z 2 , z 3 ) + 2P 3 (z 1 , z 2 , z 3 ) = — + -G 13 - ^G 123 

= ± + i G ( Zl3 )_iJ^ + ^ + £M\ (3 . 16) 

12 6 6 { z 13 z 23 z 12 z 32 z 21 z 31 J 

Expanding G-function as 

G(z) = -2(1 + z 2 ) (C(l) - z 2 ( a (3) + •••}, (3.17) 



14 



and substituting it into (13.161) . one can take the homogeneous limit z 1 , z 2 , z 3 — > 0, 

lim (S^SI) ( Zl , z 2 , z 3 ) = i - ic.(l) - (C.(l) - C«(3)) 

= ^-|ln2 + ?C(3), (3.18) 

which reproduces the expression (11.31) 

4 Inhomogeneous correlation functions for n = 4 

Next we consider the inhomogeneous correlation functions for n — 4, Pele£eie£{zi, z 2 , z 3 , z 4 ). 
The EFP Pi(zi, z 2 , z 3 , z 4 ) = P++++ (zi, z 2 , z 3 , z 4 ) was calculated in [H41. 

*4\Z\, z 2 , z 3 , , 2 2) ^3? Z4) (7(212) + 4l 1(^1, Z 2 , Z i )G(Zi 3 ) 

5 

, Z4, Z 3 , Z 2 )G(Zn) + A 4 1(23, Z 2 , Zi, Z<i)G(z 23 ) 
l{Z4, Z 2 , Z 3 , Zi )G(z 2A ) + Ay l{Z4, Z 3 , Z 2 , Zi 

+ A^(zi, z 2 , z 3 , Z4)G(zi 2 )G(z 34: ) + A^ 2 (zi, z 3 , z 2 , z 4 )G(zi 3 )G(z 2 4 : ) 
, z&, z 3 , z 2 )G(zu)G(z 23 ) , 



with 



a(0)( \ Qf}( Z ^ Z 2: z 3,Z4) 

A\J(zi,z 2 ,z 3 ,z 4 ) = , (1 = 1,2), 



Zl 3 Zi4Z 23 Z 2 4 

(o) / n _ 1 / / , w ^ , 2 / 2 3 \ / 2 3 \ 3 



Q\i(zi, Z 2 , Z 3 , Z 4 ) = — j(zi 3 Z 2 4 - 1)(214Z23 ~ 1) + j7 Ul2 + g ) ( Z 34 + 2 j ., 
Q 4 °J(z 1; 22, Z 3 , «4) = 2Q$(*l, 22, Z 3 , * 4 ) - ^(4 + 4) (z\\ + 1). (4.1) 

Let us think of the correlation functions in the sector e x + e 2 + e 3 + e 4 = e' x + e 2 + e 3 + e 4 = 2. 
There are 16 correlations in this sector. They, however, are connected via the intertwining relations 
similarly in the case of n = 3. In fact, if one can calculate a single correlation, for example, 
P+++^(zi, z 2 , z 3 , z 4 ), other correlation functions are obtained with the help of the intertwining 
relations and the reduction relation 

P-++- + ( z i, z 2 , z 3 , z 4 ) = P 3 (z 2 , z 3 , Z4) — PAzi, z 2 , z 3 , Zi). (4.2) 
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Therefore we can invoke the same strategy as was done for the calculation of P++*(zi, z 2 , z 3 ) in 
the previous section. Namely, we reduce the first principle relations to those for a single functions 
P+++^(zi, Z2, Z3, 24) and solve them. Below we list the equations corresponding to (I3.5l) - d3.9l) for 

n = 3. 

• symmetry relations 

P+++- (zi, Z2, z 3 , Z4) = P + + + ^(zi, z 3 , z 2 , Z4) = P + ++^(z4, z 3 , z 2 , z\) (4.3) 

• transposition relation 

P+++- (— Zi, —z 2 , —z 3 , —Z4) 

_ (Zl2 + i)(zi3 + Z)(Z43 + i)(*43 + l) v 

— " +++ _ ^lj z 2 , z 3 , Z A ) 

Z\ 2 Z\ 3 Z A2 Za 3 

^12 + i)(z 14 + i)(z 32 +i) p _ +++r x 

^12^14^32^43 

i{z l3 + i)(z u + i){z 23 + i) „_+++, x 

[Zl, Z 3 , Za, z 2 ) 

z\ 3 z\az 23 z A2 
i{z43 + i){z4i +i)(z23 +i) 



P + ++-(z 2 , 2l, #3, £4) 



ZA3ZA1Z23Z12 

(^3, ^1, ^2, 2:4) 



i(z A2 + i)(z 41 + i)(z 32 + i) n _ + 



Z A2 Za\Z 32 Z\ 3 

1 + i (zi 2 + z A3 ) + z 12 z 3A + z 13 z 2A +++ 

H "+++_ \ Z 3i Z l) z ii z ±) 

Z\ 2 Z\ 3 Z 2A Z 3 a 

+ ifa 2 + i)fa + i) P3(2i22| 2s) + ^ + 0(-- u + ftfa| , 3 _ Zi) + JE±L PM , 22 , 24 ) 

^42^43^41 Z\ 2 Z\ 3 Z\4 Z 32 Z A3 Z 3 \ 

+ ^ 23 + Z P 3 (z 2 , 23,24) (4.4) 

^23^42^21 

first recurrent relations 

lim Prit"j"(^i, ^2, ^3, 24) = 0, 

lim P~+++(2i, z 2 , z 3 , 24) = -P+^(z!, z 3 , z A ) (4.5) 
second recurrent relation 

P+++.t(z + i, z, z 3 , 24) = 0, 

P+U-( z ~ h z , -23, z 4 ) = -P-++{z, z 3 , z A ) (4.6) 
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identity relation 



-242 + i z 43 + i 



P+++-(zi, z 2 , z 3 , z A ) H P+++_(zi, z 3 , z±, z 2 ) 



Z42 Z43 z 23 z 2 4 

+ Z32 + % ZU + i P+++- (zi, z 2 , z±, z 3 ) = 2P 4 (zi, z 2 , z 3 , Zi) - P 3 (z 2 , z 3 , Z4) (4.7) 

^32 ^34 

• translational invariance 

P+++-(Zl +X,Z 2 + X, Z 3 + X, Z 4 + X) = P^+++(2i, Z 2 , Z 3 , Z4) (4.8) 

According to the ansatz (12.141) . (12.231) . P+++j~(,zi, z 2 , z 3 , z 4 ) may be written in the form 

P+++-{zi, z 2 , z 3 , Z4) = — + A^\(z\, z 2 , z 3 , Z4)G(zi 2 ) + A^(zi, z 3 , z 2 , Zi)G(zi 3 ) 
+ A^ 2 (zi, Z4, z 2 , z 3 )G(zu) + ^4^3(^3, z 2 , z\, Zi)G(z 23 ) 

+ 4 1(^4 , Z 2 , Z 3 , Zi)G{z 2 4) + A 4 i(z4, Z 3 , Z 2 , Z\ )G(z u ) 

, z 2 , z 3 , Z4)G(zi 2 )G(z 3 4) + A A 4 (zi, z 3 , z 2 , z i )G{zi 3 )G{z 2 4) 
+ 44,5( 2 l5 ^4, z 3 , z 2 )G(zi4)G(z 23 ), 

AWf \ QaJ( Z ^ Z ^ Z 3,Z4) 

A\l(z l} z 2} z 3 ,z 4 ) = — '■ , (z = l,...,5), (4.9) 

^13^14^23^24 



where we have taken the symmetry relations (14.31) into account. We can assume the polynomials 

Qij (zi, z 2 , z 3 , Z4) as 

Q\%^z 2 ,z 3 ,Z4) = C±AhJ2j 3 ,k)zi 1 zi 2 4 3 z{\ (4.10) 

0<Ji .32 J3J4< 2 
il+J2+J3+J4^ 4 

Substituting (14.91) and (14.101 ) into (14.41) ( 14.81) . we can extract the equations for coefficients C^ji, j 2 , j 3 , j 4 ) 

and solve them. Actually we have made a Mathematica program to perform the calculation. The 
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result is 



Qi}( z l} z 2i z 3i z i) 

19 7 2 1 2 ^ 2 _ 1 2 1 2 ^ 2 _) 2 2 , ^ 2 2, _J 2 2 

" 90 + 8(f 12 12^ 3 144 ^ 14 36 * 23 48* 24 + 30 ^ 144 V24 + 144 V23 + 720 Zl2Z34 

, Y 1 1 1 2 1 2 1 2 ,5 2 1 2 A 

V 9 Zl3 ~~ 36 Z2A 72 212234 _ 18 212234 ~ 36 2l3224 72 213224 _ 24 214223 J > 

Q^^l) Z 2, ^3) ^4) 

19 59 2 ^~ 2 2 ^~ 2 ^~ 2 ^2 ^~ 2 2 ^~ 2 2 17 2 2 

= 90 + 360* 12 ~ 12 Zl3 ~ 12 * 23 ~ 12^ 14 " l^ 24 + 45^ + 18 ^ 24 + 18 V23 " 360 V34 

1 1 2 1 2 1 2 1 2 \ 

V 9 Zl3 9 224 ~~ 18 213224 ~ lg Zl3Z2i ~ J^ Z U Z 23 - Y^ Z U Z 23 I ) 
QiJ^l) z 2, ^3) 24) 

19 29 2 ~^ 2 2 ~^2 2 17 2 ^~ 2 2 ^~ 2 2 ^~ 2 2 

= 90 + 360 212 " is^ 13 " 48 ^ 23 " 48^ " 48* 24 + 360* 34 + 72* 13 * 24 + 72 * 14 * 23 " 180* 12 * 34 

+ z 213 

V 36 36 



. 1 2 1 1 2,-^2 I 21 

■^24 + -^13^24 + ^^13^24 + ^^14^23 + J^ Z U z 2 3 J ' 



Qa\( z 1, z 2i z 3i z a) 



7 1 2 J_ 2 _J_2 ]_ 2 J_ 2 J_2 _J 2 2 _ J_ 2 2 , ;j_ 2 2 

~ 45 + 15 2 " 12 24 Zn 72 Zu 72^ 24^ 24 + 15 ^ + 360* 12 * 34 72 213224 + 72 * 14 * 23 

L Y1 1 1 2 1 2 1 2 1 2 

\ 18 213 ~~ 18 224 36 212234 _ 36 212234 _ 18 213224 18 213224 

^4,5 (^1 5 z 3i z i) 

"^" 2 ~^2 ^~ 2 ~^2 ^~ 2 17 2 ^~ 2 2 ^~ 2 2 ^~ 2 2 

= 45 + 180 212 ~ 24 %3 ~ 24 214 " 24 * 23 " 24 224 + lSO^ 34 " 90 ^ 34 + 36 * 13 * 24 + 36* 14 * 23 

, Y 1 1 12 1 2 1 2 1 2^ „ m 

+ \l8 Zl3 + 18 ZM ~ 36* 13 * 24 ~ 36^ 24 ~ 36 V23 ~ 36* 14%J ) ' (4,11) 

Unfortunately we could not find a further compact way to express the results above. Other 
correlation functions in this sector are calculated by the intertwining relations. As it is not pos- 
sible to present all the explicit results here, we instead show the formal expressions in Table 
2. Here C(zi, z 2 , z 3 , Z4) and Z(z±, z 2 , z 3 , 24) are respectively the real and the imaginary part of 
the correlation function P7^++*(zi, z 2 , z 3 , z 4 ) (14 .lib . Similarly, A(z\, z 2) z 3) z A ), B(zi, z 2 , z 3 , z 4 ), 
D(zi,z 2 ,z 3 ,z 4 ), E(z 1 ,z 2 ,z 3 ,z A ) and X(z u z 2 , z 3 , z 4 ), Y(z u z 2 , z 3 , z A ) are respectively the real 
and the imaginary parts of other correlation functions. Using the intertwining relations, they are 
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p e l e 2 e 3 e 4 
- r eie2e3£4 


- + + + 


+ - + + 


+ + -+ 


+ + +- 


- + + + 


Pz\Z 2 ,Z 3 ,Z 4 ) 
-P 4 (z 1 ,Z 2 ,Z 3 ,Z 4 ) 


A ( \ 

A(z 1 ,Z 2 ,Z 3 ,Z 4 ) 

-iX(z 1 ,z 2 ,z 3 ,z 4 ) 


B{z 1 ,z 2 ,z 3 ,z 4 ) 
-iY(z 1 ,z 2 ,z 3 ,z 4 ) 


C{Z 1 ,Z 2 ,Z 3 ,Z 4 ) 

-iZ(z 1 ,z 2 ,z 3 ,z 4 ) 


+ - + + 


A I \ 

A(z 1 ,Z 2 ,Z 3 ,Z 4 ) 

+iX(z!,z 2 ,z 3 ,z 4 ) 


D(zi, z 2 , z 3 , z 4 ) 


E[Zl,Z2,Z 3 ,Z 4 ) 

(^123 + ^234) 


B(z 4 ,z 3 ,z 2 ,zi) 
-iY(z 4 ,z 3 ,z 2 ,zi) 


+ + - + 


B{Z 1 ,Z 2 ,Z 3 ,Z 4 ) 

+iY(zi,z 2 ,z 3 ,z 4 ) 


TP/ \ 

E(z 1 ,z 2 ,z 3 ,z 4 ) 

~~k 4 (^123 + QiW) 


D(z 4 ,Z 3 ,Z 2 ,Z 1 ) 


A ( \ 

A(z 4 , z 3 , z 2 , zx) 
-iX(z 4 ,z 3 ,z 2 ,Zi) 


+ + +- 


C(Z 1 ,Z2,Z3,Z4) 
+iZ(Zl,Z 2 ,Z 3 ,Z 4 ) 


B{z 4 ,Z 3 ,Z2,Z\) 

+iY(z 4 ,z 3 ,z 2 ,zi) 


A(z 4 , Z 3 , Z 2 , Zi) 

+iX(z 4 ,z 3 ,z 2 ,zi) 


Ps{zi,Z2,Z 3 ) 

-P 4 (zi,z 2 ,z 3 ,z 4 ) 



Table 2: Inhomogeneous correlation functions for n = 4 in the sector ei + e 2 + e 3 + e 4 = 2 

successively calculated from C(z 1: z 2 , z 3 , z 4 ) and Z(zx, z 2 , z 3 , z 4 ), 

B(zi, z 2 , z 3 , z 4 ) = C(zi, Z2, z 4 , z 3 ) H {Z (zi, z 2 , z 3 , z 4 ) — Z(zi, z 2 , z 4 , z 3 )} , 

z 34 

Y{Z\, z 2 , z 3 , z 4 ) = Z{zi, z 2 , z 4 , z 3 ) H {C(zi, z 2 , z 4 , z 3 ) — C(zi, z 2 , z 3 , z 4 )} , 

^34 

A{z u z 2 , z 3 , z 4 ) = B{z u z 3 , z 2 , z 4 ) H {Y{zi, z 2 , z 3 , z 4 ) - Y{z x , z 3 , z 2 , z 4 )} , 

£23 

X[Z\, z 2 , z 3 , z 4 ) = Y(z\, z 3 , z 2 , z 4 ) H {B(zi, z 3 , z 2 , z 4 ) — B[Zi, z 2 , z 3 , z 4 )} , 

Z23 

D{Zi, Z 2 , Z 3 , Z 4 ) = P 3 (zi, Z 3 , Z 4 ) — P 4 (z 2 , Zi, Z 3 , Z 4 ) H {X(z\, Z 2 , Z 3 , Z 4 ) — X{Z2, Zi, z 3 , z 4 )} , 

^21 

Efa, z 2 , z 3 , z 4 ) = B(z 2 , z lt z 3 , z 4 ) H {Y{z x , z 2 , z 3 , z 4 ) - Y(z 2 , z u z 3 , z 4 )} . (4.12) 

£21 

Note that the diagonal correlation functions ( z u z 2, z 3 , z 4 ) and P++I+(zi, z 2 , z 3 , z 4 ) have 

no imaginary parts. Note also that we have found the imaginary part of P+7*+(zi, z 2 , z 3 , z 4 ) is 
simply expressed in terms of the function G\\ 3 introduced in (13.11) . These properties are reflected 
in some symmetry relations of the functions A(zi, z 2 , z 3 , z 4 ),..., Z(zi, z 2 , z 3 , z 4 ). For example we 
have 

, z 2 , z 3 , z 4 ) = A(z 2 , Zi, z 3 , z 4 ), E(zx, z 2 , z 3 , z 4 ) — E(z 1 , z 3 , z 2 , z 4 ), 

Y(z±, Z 2 , Z 3 , Z 4 ) = Y(z 2 , Zl, Z 3 , Z 4 ) + — (G^23 + ^234 — ^213 — ^34) ■ (4-13) 
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p e l e 2 e 3 £ 4 
- r eie2£3£4 


i 

-\ — h 


i i 
-t- -t- 


-r ~r 


i 

T- + 


1 1 

-r + 


1 

-r -r 


1 — h 








CJS ■y-W 


U 5U— z ij) 


r*„(l -y-W 


1 r 






^SKi-Zij) 


Lj 9\\-Zi\) 


C<,„(i -y-W 
L-lOU^i i) 


U 5K\Zi}) 


I i 

n r 






^liKlZif) 


^12\\ Z %S) 






+ h 






Cn({z t }) 


Cu({Zi}) 


Cs({zi}) 


Cs({zi}) 


+ -+- 


C,({zi}) 


C 10 ({zi}) 






C 7 ({zi}) 


C2({Zi}) 


+ + 


Ce({zi}) 


C 5 ({-Zi}) 




C 3 ({-Zi}) 


C2({-Zi}) 


Ci({zi}) 



Table 3: Inhomogeneous correlation functions forn = 4 in the sector e\ + e 2 + e 3 + e 4 = 



We also remark the following relations hold from the identity relations. 

A(zi, z 2 , z 3 , Z4) + B{zi, Zi, Z3, Z4J + C(zi, Z2, z 3 , Z4) = 2P 4 (z 1 , Z2, z 3 , Z4) — P 3 (z2, z 3 , Z4), 
A(zi, Z2, z 3 , Z4) + D{zi, z 2 , z 3 , Z4) + E(zi, z 2 , z 3 , Z4) + -8(2:4, z 3 , z 2 , z\) = P^zi, Z2, z 3 , Z4), 
X(z±, z 2 , z 3 , Z4) — Y(z4, z 3 , z 2 , zi) — — (Gi 23 + G234) = 0, 

X(zt, z 2 , z 3 , Z4) + Y(z!, z 2 , z 3 , Z4) + Z(zi, z 2 , z 3 , Z4) = 0. (4.14) 

Unfortunately we have not been able to afford a nicer compact way to describe the functions 
A(zi, z 2 , z 3 , Z4),..., Z(zi, z 2 , z 3 , Z4), from which the relations such as (14.131) . (14.141) follow natu- 
rally. 

Next let us consider the sector ei + e 2 + e 3 + e 4 = e[ + e' 2 + e' 3 + £4 = 0. There are 36 correla- 
tion functions in this sector. But according to (12.81) . we need to know only 12 correlation functions 

[z\z2i z 3 , Z4),...,Ci2(zi, Z2, z 3 , Z4) as shown in Table 3. 
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In the present case, they can be derived from those we have obtained already as follows. 



C\{Z\, Z 2 , Z 3 , Z4) — P_++(z 2 , Z 3 , Z4) — P+_++(Zi, Z 2 , Z 3 , Z4), 

11 1 

= 12 12 ^ 34 ~~ ^ 23 ~ ^ 12^ 234 _ ^( Zl > z<2 > z?n Zi 'i 

C 2 (zi, z 2, z 3, Z4) = P + ^+(z 2 , Z 3 , Z4) — P+ + ^+(zi, Z 2 , Z 3 , Z4), 

= ^ + Q^i2 - E(z l , z 2 , z 3 , z 4 ) + — (Gf 23 - G 2U ) , 
, z 2 , z 3 , Z4) — P + + + (z 2 , Z 3 , Z4) — P+ + +*(zi, Z 2 , Z 3 , Z4) 

= J^ + g^ 24 ~ g G 23 4 ~ B ( Zj i, Z 3, z 2 , Zi) - i ^Y(z 4 , z 3 , Z 2 , Zi) - j^G 23 4 }> , 
C4( Zl , Z 2 , Z 3 , Z4) = C 3 (-Z 4 , ~Z 3 , -Z 2 , - Zl ) 

= T2 + \° 13 ~ ^ Gl23 ~ B(yZl ' Z2) Z3 ' Z ^ + * { Y ^ Zu Z2) Zs ' Z ^ ~ T2 G ^ 
65(21, z 2, Z 3 , Z4) = ( 1 H ] C 3 (z 2 , Zi, z 3 , z 4 ) C 3 (zi, Z 2 , Z 3 , Z4), 

V Z 12 J z 12 

5 

Cq(zi, z 2 , z 3 , Z4) = Pa(zi, z 2 , z 3 , z 4 ) — Cj(zi, z 2 , z 3 , Z4), 

3=1 



Cj(z\, Z 2 , Z 3 , Z4) — -P_+_ ( z i, Z 2 , Z 3 ) — P_+ (Zi, Z 2 , Z 3 , Z4) 

= P+-+ ( z l, z 2, Z 3 ) — P+-++(zi, Z 2 , Z 3 , Z4) 

11 1 

= 12 12 ^ 13 ~~ ^ 23 ~~ ^ 12 ^ ~~ l^^ 123 ~~ ^( Zu Z2) Zs > ^ 4 )' 

, z 2, z 3, z a) — P+ + -( z 2, z 3, z i) — P+++-( z l, z 2, z 3, z i) 

= ^+ g G 34 - A(Z4, Z 3 , Z 2 , Zi) - i ^X(Z4, Z 3 , Z 2 , Zi) + ^ G 234 } , 
C 9 (zi, Z 2 , Z 3 , Z4) = C 8 (-Z4, -z 3 , -z 2 , -Zi) 

= ^+ g^ 12 ~~ M Z 1, Z 2, z 3, z a) + i ^( z l, z 2, z 3, z i) + J^^123 J ' 

1 

Cio(zi, z 2 , z 3 , Z4) = Cq(zi, z 3 , z 2 , Z4) H {C§(zi, z 3 , z 2 , Z4) — C$(zi, z 2 , z 3 , Z4)} , 

Z32 

CiM 

, z 2, z 3, z 4) —++ \ Z 1, Z 2, z 3) ~ PI +++ ( z l, z 2, z 3, z i) 

= —- — — (Gi 2 + G13 + G 2 4 + G34) + — (G123 + G234) + Pa( z 1, z 2, z 3, z a), 
12 12 

Ci 2 (zi, z 2 , z 3 , z 4 ) = C w (zi, z 2 , Z4, z 3 ) H {C 9 (zi, z 2 , z 4 , z 3 ) - C 9 (zi, z 2 , z 3 , z 4 )} . (4.15) 

£43 
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Among the correlation functions above, let us pick up a diagonal one P ++ (21, z 2 , z 3 , z 4 ), 

since it has a relatively simple explicit form. 

P--++ (Zl, *2, Z3, Za) = ^ + Af^, Z 2 , 2 3 , Z 4 )G(z 12 ) + |aJ 2 J(zi, Z 3 , Z 2 , 2 4 ) ~ j G(z 13 ) 

+ |^4 2 i(^i> z 4 , z 3 , z 2 ) - j G(z 14 ) + |j4$(z 3 ,Z2, 21,24) - j £(223) 

+ I A? }(z 4 , 2 2 , Z 3 , Zi) - T^r 1 G(2 24 ) + A?J(Z 4 , Z 3 , Z 2 , 2i)G(z 34 ) 



12 , 

+ A^(zi, 2 2 , z 3 , z A )G(zi 2 )G(z 3A ) + A^ 2 (zi, z 3 , z 2 , z 4 )G(z 13 )G(z 24 ) 
+ A 4 l(zi, z 4 , z 3 , z 2 )G(z 14 )G(z 23 ), 



with 



a{2)( \ Q ( S{ZI,Z 2 ,Z 3 ,Z 4 ) 

A\l{z u z 2 ,z 3 ,z 4 ) = , (2 = 1,2), 

2l32l 4 ^2322 4 



3\ 5 
2 J ~ 2 



Q 4 2 2(2l, 2 2 , Z 3 , 2 4 ) = |(2l 3 2 24 + l)(z 14 Z 23 + 1) + ^ ^ 2 + ^ 
QS(2i,Z 2 ,2 3 ,Z 4 ) = 2Q^(zi, 22,^3,^4) - ^j(*12 - !)(^34 + !)■ ( 4 - 16 ) 



The correlation functions in the remaining sector ei + e 2 + 63 + e 4 = e'x + e' 2 + e 3 + €4 = — 2 are 
obtained from those in e x + e 2 + e 3 + e 4 = e' : + e 2 + e 3 + e 4 = 2 by the negation relation 

piltf (21, 2 2 , z 3 , Z4 ) = p::i;::^;:::( Zl , Z2 , Z3 , Z4 ). had 

In this way we can get all the inhomogeneous correlation functions for n = 4. 

Now we shortly consider the diagonal elements of the correlation functions P^tztt { z i, z 2 , z 3 , z 4 ). 
We claim they can be expressed only in terms of the EFP's, P 2 (z 1: z 2 ), P 3 (2i, z 2 , z 3 ), P 4 (z 1: z 2 , z 3 , z 4 ) 
and PIT++ (zi, z 2 , z 3 , z 4 ). For example, we have 

P-+++ (2l, 2 2 , Z 3 , Z 4 ) = P 3 (z 2 , Z 3 , Z 4 ) — P 4 (Zi, Z 2 , Z 3 , Z 4 ), 

(2l, 22, 2 3 , Z 4 ) = P2(2 3 , 2 4 ) — P 3 (^2, 2 3 , Z 4 ) — P ++ (21, ^2, 2 3 , Z 4 ) , 

1 

2 



^+_+_ (2l, 2 2 , Z 3 , Z 4 ) = - - P 2 (2l, 2 2 ) - P 2 {z 2 , Z 3 ) - P 2 {z 3 , Z 4 ) 



+ P 3 (2l, 2 2 , Z 3 ) + P 3 (2 2 , Z 3 , Z 4 ) + P__++ (21, 2 2 , Z 3 , 2 4 ), 
P+--+ (2l, 2 2 , 2 3 , Z 4 ) = P 2 (2 2 , 2 3 ) - P 3 (2i, Z 2 , 2 3 ) - P 3 (2 2 , Z 3 , Z 4 ) + P 4 (2i, Z 2 , Z 3 , Z 4 ). (4.18) 
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Then, especially, the inhomogeneous generalization of the third-neighbor correlator (S^S^) is rep- 
resented 



{Si SI) \Z\, Z2,Zs, Zi) 

= ^|^++++ (^1) z 2, Z3, Zi) — P_+++(zi, Z 2 ,Z 3 , Z4) + P+_++(zi, Z2, Z3, Z4) — P — ++(zi, Z2,Z 3 , Z4) 

+ P++-+[zi, Z2, z$, Z4) — P_ + _ + {z\, Z2, Z3, Z4) + P+ + (zi, Z2, Z3, Z4) — P +(zi, Z 2 , Z3, Z4) 

— P+++- ( z l, Z2, Z3, Z4) + P_+ + _ {Z\, Z 2 , Z3, Z4) — P+_+_ (Zi, Z 2 , Z 3 , Z4) + P — + _ (Zi, Z2, Z3, Z4) 

— P++ (zi, z 2 , Z3, Z4) + P_+ (zi, z 2 , z 3 , Z4) — P^ (zi, z 2 , Z3, Z4) + P (zi, Z2, Z3, Z4) I 

= -| + z 2 ) + P-iiz-i, Z3) + P 2 (z 3 , z 4 ) - 2 {P 3 (zi, Z 2 , Z3) + P 3 (z 2 , Z3, Z^} 

+ 2 {P 4 { Zl , Z 2 , Z3, Zi) — P — ++(z\, Z 2 , Z3, Z4)} 

— — + -Gu + - (G123 + C234) + A^l(zi, Z2, Z3, Zi)G\2 + Afl(zi, Z3, z 2 , ^4)^13 
12 o 

, Zi, Z3, Z 2 )Gn + A 4 i(z3, Z 2 , Zi, Zi)G23 + A 4 i(Z4, Z 2 , Z3, Z\)G2i 

, Z3, Z2, zi)G 3 i + A 42 { z i, z 2 , Z3, Zi)GuG3i + A 4 2 { z i, Z3, z 2 , £4)^13^24 

, Zi, Z3, Z2)GuG23, 

where 

|( 3 )^ „ „ ,\_oJ/|(»)/, „ „ „\ 4(3), 



) z 2, Z3, Zi) = 2 {A%( Zl , z 2 , Z3, Zi) — A 42 (zi, z 2 , Z3, Zi) 

I ( I ! I.) , 

9 \Z\3Z\iZ23Z2i Z\iZ23 Z\3Z2i) 

) Z 2 , Z3, Zi) = 2{A{% 1 

) z 2, Z3, Zi)-A%[{ Zl , Z 2 , Z3, Zi) 
1/5 1 1 



18 \Z13ZnZ23Z2i Z\iZ23 Z 13 Z 2 i Z\iZ2i Z13Z23 

By taking the homogeneous limit z± — > 0, we recover the result (I1.6I) . 

lim (SJS*) (21, z 2 , Z3, Zi) = ±-- ic.(l) + 2 {C.(l) - Ca(3)} 



+ { - yCa(l) + f C.(3) - yC.(5) " f Ca(l)Ca(3) - j|Ca(3) 2 + yC(l)C.(5) 



(4.19) 



fc u (l)C(3)-^(3) 2 + f< 

= ^ " 3Ca(l) + yCa(3) - f Ca(5) - yC(l)C(3) - ^Ca(3) 2 + f C(l)Ca(5) 

= - - 3 In 2 + ^C(3) - -^C(5) - y m 2 • C(3) - §C(3) 2 + -| In 2 • C(5). (4.20) 
Other correlation functions in (I1.6I) are reproduced in a similar way. 
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5 Correlation functions in the case n = 



5 



Following the scheme in the previous sections we could further obtain all the correlation functions 
for the XXX chain on five lattice sites. In fact we have calculated two inhomogeneous correla- 
tion functions z 2 , z 3) z 4 , z 5 ) and P+++^*(zi, z 2 , z 3 , z 4 , z 5 ) from the first principle re- 
lations. Other correlation functions are derived by the intertwining relations, etc., as before. The 
calculations have been performed by the heavy use of Mathematica and the obtained results are 
too complicated to be described thoroughly in this paper. Therefore here we shall omit to describe 
them. We, however, give explicit form of the inhomogeneous generalization of the fourth-neighbor 
correlation function (Sf )(z±, z 2 , z 3: Z4, z 5 ). 

(SiS^)(Zl, Z 2 , Z 3 , Z4, Z 5 ) = - ^2 e l e 5^e?£2 2 e3 i e?e5 5 ( Zl ' ^ 2 > Z ^ Z $) 

= — + ^.5,i(^i) z 2 , z 3 , Z4, zs)Gi 2 + ^1(21, z 3 , z 2 , Z4, zn)Gi 3 + A^i(zi, Z4, z 3 , z 2 , za)Gi4 

+ ^-5,1 ( z 5, z 2i z 3i z 4, z l)G 2 s + A 5jl (z$, Z3, Z 2 , Z4, ^1)^35 + A 5> i{Z5, Z4, Z3, Z 2 , ^i)G45 

+ As j2 (z 2 , z 3 , zi, Z4, z^)G 23 + A^ j2 (z 2 , Z4, zi, Z3, z^)G 2 4 + A^ j2 (z 3 , Z4, z±, z 2 , z^)G 3 4 

+ ^5,3(^1) z 5, z 2, z 3i ^4)^15 + As t 4{Zi, Z 2 , Z3, Z4, Z$)Gi 2 G34 + A^^Zi, Z3, Z 2 , Z4, ^5)^13(^24 
+ ^5,4(^1? z 4, z 2i z 3i Z ^)G\4G 23 + A 5i4 (z5, Z 2 , Z3, Z4, Zi)G 2 $G34 + ^4(25, Z3, Z 2 , Z4, Zi)G 3 ^G 2 4 
+ ^5,4(^5; z 4, z 2, z 3, ^1)^45(^23 + ^.5,5(^1, z 2, z 3i z 5, z i)Gi 2 G35 + A^^(zi, Z3, Z 2 , Z5, £4)^13(^25 
+ ^5,4(^1, z 2i z 4, z 5, £3)^12*^45 + ^5(21, Z4, Z 2 , Z5, Z 3 )Gi4G 2 5 + A^^Zi, Z3, Z4, Z5, Z 2 )Gi 3 G45 

+ ^5,4(^15 z 4i z Zi z hi ^2)^14^35 + ^5,5(^1; z h-, z 2i z zi z 4)G\z ) G 23 + A 5i5 (zi, Z5, Z3, Z4, z 2 )Gi^G 3 4 

+ ^5,5(^1? z hi z 2-> z 4i z 3)GisG 2 4, 



where 

Qh,i{ z l, z 2, z 3i z 4, z h) 

A^i(zi, z 2 , Z3, Z4, z 5 ) = 



1 

^13^14^15^23^24^25 



1,2,3) 



Qb,i( z l, z 2, z 3, z 4, z 5) /. r x 

-, (2 = 4,5) 



, ^13^14^15^23^24^25^34^35 
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with 

$5,1 (^lj z 2, z 3, Z4, Z 5 ) = — [(z 15 + Z 2 5)(zi 3 + Z 2 a){2 — Z 13 Z U — Z14Z23) ~ 6(zi 3 Z u — 2)(z 2 3Z 24 — 2) 

+ 2i 2 (^13^23 + ^14^24 - 16) + lO^ - 2] , 

Q5,2(2l, Z 2 , Z3, Z4, Z5) =Q5,l(2l, Z 2 , 23, Z4, Z5) + —2i32; 2 3 [5z u Z 2 4 + 245(^14 + Z 2 4) — 5] , 

lo 

$5,3(2l, Z 2 , Z3, Z4, Z5) =$5,l(2l, Z 2 , Z3, Z4, Z5) + — Z15225 \&Z\3Z\4Z 2 3Z 2 4 — A(z\3Z 2 4 + Z^Z^) — 3z| 4 

$5,4(21, Z 2 , Z3, Z4, Z 5 ) —— [2(^13^24 + 214^23)^15^25 — 213^23^24245 _ 214^23^24235 — 213^14^23^45 

lo 

1 2 2 1 a( 2 1 2 \ d / 2 , 2 1 2 1 2 \ 

— ^13^142:24235 + 2i 2 2 3 4 + t>(2 12 + 2 34 J — ~ ^3 + Z u + £ 23 + £ 24 J 

- (2l 2 5 + 2 2 2 5 + 4 + 4) + 12] , 

$5,5(21, z 2 , Z3, Z4, Z5) =$5,4(21, z 2 , Z3, Z4, Z5) + — (2 — z 13 z 24 — 214^23)^15^25235245- (5.21) 

By taking the homogeneous limit, we obtain the fourth-neighbor correlation function for the ho- 
mogeneous Heisenberg model as follows. 

(SjSj +4 ) = lim (S'S§) (zi, z 2 , z 3 , Z4, z 5 ) 

= Y 2 - f c ° (1) + 1T U3) - 72C « (1) ' c ° (3) - T c " (3)2 - ir c " (5) 

+ ^ui)-U5)-^c.(3).u5)-fc^ + fa7) 

- ^fUi) ■ CP) + 280C„(3) • C(7) 

^-f'" 2+1 f« 3 »- 541 " 2 '« 3 »-^( 3 ) 2 -^( 5 ) 
+ lip in2 .C(5)-f C (3).C(5)-^ + ^C(T) 

1715 6615 
— In 2 • C(7) + — — C(3) • C(7) = 0.034652776982 • • • . (5.22) 

This is one of the main new results of this paper. Furthermore we have obtained all the other 
correlation functions for n — 5. We present below the homogeneous limit of independent ones. 

(S*S* +1 S !+2 S> +t ) = 4 + ^2 - «Z«3) + 1^2. C (3) + ^«3)' + f| C(5) 

_^ ln2 . ((5) + ^ ( 3). ( (5) + f ( (^-|5 ( (7) 

1029 441 
+ — In 2 • C(7) - — C(3) • C(7) = -0.0098924350847 • • • , (5.23) 
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1 1 301 Q 1fl7Q Q75 

(S 1SU S* +2 S- +l ) - Wa --^2 + - C (3) - - in 2 . ((3) - — { (S)' - -«5) 

365, n . . 185 A/oN . . 1375 A/rx2 735 

+ - 2 ■ C(5) " I28« 3) ' C(5) " T2¥ C(5) + I28 C(7) 

1323 4851 
- — In 2 ■ ((7) + ^gg-C(3) • C(7) = 0.0027889733995 • • • , (5.24) 



+ f.n2.C, 5 ,->,.C( 5 »-^ + Ic(7, 

i Qo^ 4851 

- — In 2 ■ ((7) + "25gC(3) • C(7) = -0.00505666089481 • • • , (5.25) 

<^J +1 ^> = ^ + iln2- ^ C (3) + f in 2 ■ C(3) + f«3)' + |f C(5) 

-^-•C( 5 ) + >-C( 5 ) + ^C(5,-^C(7, 

833 5733 
+ — In 2 • C(7) - "25g-C(3) • C(7) = 0.01857662093837 • • • , (5.26) 

(S^SM = ^ - A ln2 + ^ C (3) - f ln2 ■ C (3) - ^CO) 2 - ^C(5) 

+ 11^ i n 2 . ((5) - — C(3) • C(5) - — C(5) 2 + — C(7) 
12 64 w w 64 256 

- ^ In 2 • C(7) + ^C(3) • C(7) = -0.00108936897291 • • • , (5.27) 

<^ +1 ^ 4 > = ^ - i in 2 + ^ C (3) - f in 2 ■ C (3) - ^ C (3) 2 - ^C(5) 

1135 , « 165 /M , A/r , 1375 A/rN2 3577 

+ — In 2 ■ C (5) - -C(3) ■ C(5) - -^CW 2 + ^(7) 

- 51! 1„ 2 • ((7) + ^C(3) • C(7) = 0.001573361370145 • • • . (5.28) 

As a confirmation of our results, we have evaluated the numerical values of the correlation func- 
tions for the finite system size N up to 32 (Table 4). We have further applied an extrapolation to 
this data such as + Ci/N 2 + c 2 /iV 4 + c 3 /iV 6 , and obtained the estimated values in the limit 
iV — > oo. One can clearly observe our analytical results coincide with the extrapolated values with 
extremely high accuracy (more than 4 digits). 
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Table 4: Numerical values of correlation functions for finite systems 





A/ =26 


A T O O 

JM =28 


A T OA 

JM =30 


A T O O 

iv=32 


/AT \ 

Coo (-/V = (X)) 




0.0357233 


0. 0355713 


0.0354497 


0.0353508 


0.0346535 


1 L >XJ L tXj L *J L As \ 

^ + 1^+2^+4/ 


-0.0099654 


-0.0099553 


-0.0099471 


-0.0099405 


-0.0098925 




0.0028572 


0.0028477 


0.0028400 


0.0028337 


0.0027890 


/ ca; C2 \ 
\ j J+l J+2 j+4/ 


-0.0052334 


-0.0052083 


-0.0051882 


-0.0051719 


-0.0050564 


/ ca; ca; Qz Qz \ 
\°J °J+l j+3°j+4/ 


0.0185725 


0.0185731 


0.0185737 


0.0185741 


0.0185766 


/ Cx Cz Cx Cz \ 

\ D i D i+i D i+3 D i+4/ 


-0.0011377 


-0.0011309 


-0.0011254 


-0.0011209 


-0.0010893 


/ cx cz cz cx \ 
\ D J D J+l D j+3°j+4/ 


0.0016526 


0.0016413 


0.0016323 


0.0016250 


0.0015732 



We make several comments on our results. First, the two other diagonal correlation functions 
are given as a sum of non-diagonal ones similarly in the case of four lattice sites, 

I CZ cz cz cz \ / QX QX cz CZ \ | / QX QZ QX cz \ i / QX cz cz QX \ 

^+1^+2^+4/ - \ D i ^+1^+2^+4/ + \ D i ^'+1^+2^+4/ + \ D i ^+1^+2^+4/ 

= L _ 1 k 2 + 1!! C(3 ) _ 15 h, 2 . C(3 ) - ^ C (3) 2 - ™«5) 
80 12 40 w 2 sw 80 v ; 32 sw 

+ - in 2 ' C(5) - ^C(3) ■ C(5) - ^CW 2 + ^((7) 

4263 3969 
- — In 2 • C(7) + -r^-C(3) • C(7) = -0.0121601225799 • • • , (5.29) 



/ cz cz cz cz \ / cx cx cz cz \ I / cx cz cx cz \ i / cx cz cz cx \ 

Wj^j+l^j+^j+l/ - \ D j D i+l D j+3 D i+4/ + \ D j D i+l D j+3 D i+4/ + ^j+l^j+3^j+A/ 

1 2, n 1307 A/oN 77, n A/oN 1541 . , 2 2335 A/rN 

= n 2 + C(3) n 2 • C(3) C(3) C(5) 

80 3 240 SV ; 6 SV ; 80 SV ; 96 SV ; 

+ f.n,C(5,-fc(3,.C(5,-^ + ^C(r, 

- ^ In 2 • C(7) + ^^C(3) • C(7) = 0.0190606133356 • • • . (5.30) 
16 256 
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The expression for P(5) (11.51) is reproduced as 

I / QZ QZ QZ QZ \ | / Q 2 G 2 G 2 G 2 \ I / G 2 C z C z G 2 \ 

-t \DjDj +1 Dj +2 dj + 3/ -h ^Oj- dj +1 dj +3 dj +4 / -f- ^ \°j j+l°i+3 i+4/ 

1 10, 281 , , 45, >/n . 489, x2 6775 A/ . 

= n2 + C 3 n2 ■ ( 3 ( 3 2 C(5) 

6 3 24 sv ; 2 sv ; 16 SV ; 192 sv ; 

1225 w , 425 w , w , 12125 w , 9 6223 w , 
+ — 111 2 • C(5) - ^ C(3) • C(5) - ^6- C(5) + ^ C{7) 

_ in 2 • C(7) + ^^C(3) • C(7) = 2.0117259 • • • x lO" 6 . 

64 sv ' 512 sv ' sy ' 

In fact, arbitrary correlation functions on five lattice sites are similarly calculated. 
Finally we notice that the "next-nearest chiral correlator" simplifies drastically as 



((Sj x Sj + i) • (Sj + 3 x Sj+4)) — 6 ((SjSj +1 Sj +3 Sj +4 ^ (SjSj +1 Sj +3 Sj +4 j) 

= -0.015976382058- ■• , (5.31) 



as well as a similar correlator 



{(Sj x Sj + i) • (Sj +2 x Sj +4 )) — 6 {(SjSj +1 Sj +2 Sj +4 ) — (SJ Sj +1 Sj +2 Sj +4 )) 

47 7 Q 17^ 9f> 

= In 2 - -C(3) + i In 2 • C(3) + ^C(3) 2 + -^C(5) - -|m2 • C(5) 

= 0.04707380576617 • • • . (5.32) 



Recall that a similar simplification occurs in the case of "nearest chiral correlator" [34| 



((Sj x S j+ i) ■ (Sj+2 x S j+3 )) - 6 ((SjSj +1 Sj +2 Sj +3 ) - (SjSj +1 Sj +2 Sj +3 )) 

= iln2-?C(3). (5.33) 

6 Summary 

In this paper, we calculate correlation functions of the Heisenberg chain without magnetic field in 
the anti-ferromagnetic ground state using several fundamental functional relations followed from 
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the quantum Knizhnik-Zamolodchikov equations. This is the generalization of the method pro- 
posed for the emptiness formation probability. In fact we have calculated all the correlation func- 
tions on five lattice sites. Especially we could get the analytic formula for the fourth-neighbor 
correlation (SjSj +4 ). In principle we can continue to calculate (SjSj + k) k>5 > which will be re- 
ported in the further publications. 
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